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ABSTRACT 

The BATSE catalogue is searched for evidence of spindown of black holes or proto-neutron 
stars (PNS) by extracting normalized light curves (nLC). The nLC are obtained by matched 
filtering, to suppress intermediate time scales such as due to the shock break-out of GRB jets 
through a remnant stellar envelope. We find consistency within a few percent of the nLC and 
the model template for spindown of an initially extremal black hole against high-density matter 
at the ISCO. The large BATSE size enables a study of the nLC as a function of durations 
Tgo- The resulting x 2 re d is within a 2.35cr confidence interval for durations T90 > 20 s, which 
compares favorably with the alternative of spindown against matter further out and spindown 
of a PNS, whose \ 2 fits are, respectively, outside the 4a and 12cr confidence intervals. We 
attribute spindown against matter at the ISCO to cooling by gravitational-wave emission from 
non-axisymmetric instabilities in the inner disk or torus as the result of a Hopf bifurcation in 
response to energetic input from the central black hole. This identification gives an attractive 
outlook for chirps in quasi-periodic gravitational waves lasting tens of seconds of interest to 
LIGO, Virgo and the LCGT 



1. Introduction 



The nature of the inner engine of cosmological gamma-r ay bursts has been an enigma ever since their 
serendipitous discovery by the Vela (IKlebesadel et alJll973f) and Konus satellites, and even more so by 
their cosmological distribution and bi-modal distributio n in durations established by the Burst and Transient 
Source Experiment (BATSE) (|Kouveliotou et al. 1993 ). They are no w known to take place about once per 
minute throughout the Universe, tracing the earliest galaxies known ( Cucchiara et alJboilJ). At least some 
of the long GRB s are associated with supernovae of relatively heavy stars (|Paczynskill 1 99 ll : IWooslevll 1 993 : 



Pacz vnskilll998h following accurate pointing of their X-ray afterglow emission by Beppo-Sax ( Piro et al 



1995h . in a mass range that produ ces black holes or neutr on stars. For a few long GRB events, we can rule 
out an accompanying supernova (IDella Valle et all 120060 . however, and these may, instead, be associated 
with m erger events, e.g., of a neutron star with a companion black hole or another neutron star (Ivan Putten 
2008 al) . This points to a common inner engine to long GRBs with or without supernovae, whose time scale 
of duration does not appear to depend on accretion. A secular ti me scale that extends to tens of seconds may 
arise in a process of spindown of black holes or neutron sta rs dBisnovatyi-Koganlll970l : Ivan Puttenlll999i : 



Bethe et al. 2003 



van Putten 



2004: lDessartetal.l l2007. 2008). 
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Some long GRB-supernovae are hyper-energetic (iCenko et alj|2010l 120111). whose ener gy output ex- 
ceeds the maximal rotational energy of a (proto-)neutron star (PNS) (Ivan Putten et al.ll2011bh . This shows 
that PNS are not universal to all long GRBs which suggests, instead, that black holes are relevant for at least 
some of them. Distinguishing black holes and neutron stars as inner engines to long GRBs is particularly 
relevant for searches for gravitational waves associated with GRBs and core-collapse supernovae, in view of 
their rather distinct prospects for emission in gravitational radiation. Electromagnetic priors on the presence 
of black holes or neutron stars are h ereby valuable in s earches for gravitational- wave bursts by the advanced 
gravitational-wave detectors LIGO (|LIGO-Virgoll2010l) . Virgo and the Large Cryogenic Gravitational-wave 
Telescope (LCGT) kuroda et allboioh . The spin energy of rapidly rotating s tellar mass black holes ex- 
ceeds the typical energy output as observed in the electromagnetic spectrum (Frail et al.1 120011) by some 
two orders of magnitude, thus providing ample opportunity to power substantial accompanying emission 
in MeV neutrinos and gravitational waves, that latter which pointing to a sensitivity distance of about 35 
Mpc for a comple te identification of a c haracteristic quasi-periodic chirp by a time sliced matched filtering 
(TMSF) method (Ivan Putten et all 120111) . For long GRBs from mergers, an additional long duration radio 
bursts emanating from magnetic torus winds may be searched for by wide area radio survey instruments 
such as LOFAR ( van Puttenll2009 ). 



The prompt GRB emissions derive from dissipation in ultra-relativistic baryon-poor jets (BPJ) di- 
rected close to the line of sight to the observer, that hereby tracks the evolution of a long-lived inner en- 
gine by causality. The BP J may emanate from magnetic outflows from the e vent horizon of a black ho le 



Dessart et al. 



(120071 12008)). 



( Levinson & Eichler 1 1993b or from a new born proto-neutron star (PNS, e.g. 
The prompt GRB emission is hereby representative for the light curve emanating from the black hole or PNS. 
For black holes, these outflows can arise as capillary jets po wered by frame dragging, provided that the black 
hole is rotating and intermittently active ( van Puttenll2008bl ). In general terms, frame dragging introduces an 
energetic interaction between the angular momentum of particle outflow aligned with the a ngular momentum 
of the black hole. Frame drag ging has recentl y been measured by LAGEOS satellites (|Ciufolini & Pavlis 



2004 ) and Gravity Probe B jEveritt et al. 2011 ). Frame dragging induced outflows hereby establish a direct 
link between gamma-ray light curves and the evolution of black hole spin. It correctly predicts X-ray after- 
glows also from short GRBs (Ivan Putten & Ostrikeiil200ll) and points to long GRBs as a common endpoint 
of core-collapse supernovae and mergers involving rapidly spinning black holes, that may unify long GRBs 
with and without supernovae, the halo event GRB070125 and the long event GRB 050911 with no X-ray 
afterglow. 

The long duration time scales in the light curves of long GRBs hereby reflect the secular behavior of the 
GRB inner engines, after scaling out the effect of redshifts in the observed durations and different lifetimes 
of the spin of black holes or PNS resulting from unknown magnetic field strengths. For black holes, we recall 
that these are scale free objects with essentially no memory of their astronomical progenitor systems, except 
for the in itial mass and angular momentum. Black hole spindown is an alternative to spin up by accretion of 
matter (|Bardeenlll970h . where spindown is against surroundin g high density matter with accompanying low 



energy emissions that account for most of the e nergy output (|van Puttenll 19991 : Ivan Putten & Ostrikerll2001 



van Putten & Levinson 2003 



van Putten 



2008al) . The nLC is hereby sensitive to the state of the surrounding 
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matter, in particular its location relative to the Inner Most Stable Circular Orbit (ISCO). Our focus, therefore, 
will be on the asymptotic behavior in durations longer than any of the (magneto-)hydrodynamical time scales 
that are not directly relevant to the evolution of the black hole such as, in core-collapse supernovae, the time 
scales associated with black hole formation or the time scale of shock break-out of jets plowing through a 
remnant stellar envelope. In the asymptotic range of long durations, therefore, we anticipate that the nLC of 
ensemble ave rages of light curves becomes meaningful in ways similar to the Phillips light curve of Type la 
supernovae iphillipslll993h . 

The aforementioned phenomenology of long GRBs zooms in onto a specific inner engine of long 
GRBs: long GRBs with and others without supernovae point to black holes or PNS with a long duration 



time scale different from accretio n, as in thin neutrino dominated accretion disk models (iPopham et al. 



19991 : IChen & Belaborodovl I2007T ) based on the assumption of negligible output from the black hole onto 



the surrounding matter. Instead, a long duration time scale may derive from spindown, for which some 
hyper-energetic events demonstrate that a PNS is not universally viable. It is therefore natural to search for 
a prevalence of spindown of black holes in the light curves of long GRBs, since the time-averaged behavior 
in light curves should track the secular evolution of the inner engine by causality. 

For this purpose, we developed a novel time-domain analysis of the light curves of the complete BATSE 
catalogue long GRBs by application of matched filtering. The distribution of durations of long GRBs in the 
BATSE catalogue is rather broad ; which goes beyon d dilations due to their distri bution in redshift, i ndicating 
an intrinsic spread in durations (Bal azs et al.ll2003r) : confirmed by Swift, e.g.. iGehrels et al. Their 
distribution in durations, therefore, includes the diversity in intrinsic durations beyond redshift effects. The 
large number of 1491 long GRB events in the BATSE catalogue introduces an appreciable resolving power, 
that we shall exploit to extract an accurate normalized light curve (nLC) to differentiate between different 
models for spin down of black holes, with and without gravitational-wave emissions, and neutron stars. 

In §2, we review the observed high energy emissions in GRBs. §3 discusses a theory of spindown 
of black holes and PNS. §4 presents three model templates for the resulting gamma-ray light curves in 
spindown of alternatively black holes, against matter at the ISCO or further out, or PNS. In §5, we summarize 
the BATSE catalogue of light curves of long GRBs. The nLC of each model is given in §6. We interpret 
our findings in §7 for prolonged spindown of black holes differentiated for the state of matter surrounding 
the black hole, versus spindown of PNS with an outlook on long duration extragalactic bursts in radio and 
gravitational waves. Our findings are summarized in §8. 



2. Long GRBs from rotating black holes or proto-neutron stars 



The GRB light curves repr esent broad band emissions in the gamma-rays. The Fermi/GLAST event 
GRB080916C JAbdo et al.ll2009h reveals an extension of GRB spectra from 8 keV - 40 MeV of the Glast 
Burst Monitor (GBM) into the energy range 20 MeV - 300 GeV of the Large Area Telescope (LAT). 
GRB080916C is the most energetic event detected to date, of duration T90 = 66 s, z = 4.25, and with 
an isotropically equivalent energy output Ei so = 8.8 x 10 54 erg with one additional photon detected at 13 
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GeV (70 GeV in its rest frame). Here, Tqo refers to the time interval covering a 90 percentile in total photon 
count for a given burst (IKouveliotou et al.ll 19930 . Its Band parameters (Eq, a, (3) show temp oral behavior 
aroun d canonical values (500 keV, —1, —2.2) and imply a minimum Lorentz factor Fj = 890 (ITajima et al. 



l2009h . The true energy in gamma-rays satisfies 2£ 7 = f b 1 Ei so , where fa = 2/0 j denotes the beaming 
factor for a two-sided jet-like outflow with half-opening angle 9j. A number of GRB light curves display 
an achromatic break in their light curves, when the relativistic bea ming angle ~ 1 j Tj of radiation produced 



by shock fronts exceeds the half-opening angle 9j of the jet (e.g. iMeszarosI (|2006T) ). No clear measurement 



500. 



of 9j has been made for GRB080916C, but E 1 ~ 10 52 erg is expected for the typical value fa 
GRB080916C furt her reveals a delayed onset of maximal flux in the highest energy photons seen in LAT 
dTaiima et al.ll2009h . 



Prompt GRBs emissions ha ve been associated with ultra-rela t ivistic baryon poor jets emanating from 
black holes or neutron stars (e.g. Thompson ( 19941): llvletzger etal 1J201 I f), where the former is an attractive 



alternative to account for low baryon-loading (iLevinson & Eichlerj|1993"f : 



Eichler 



2011). 



If powered by the rotational energy of the black hole, they may produce essentially unlimited Eddington 
lumin osities from stellar mass black holes for the duration of black hole spin. As capillary jets (Ivan Putten 
2008bl) produced by frame dragging, they represent a direct link to the spin of the black hole. (This process 
is distinct from the Penrose process, as it operates well outside the ergo sphere.) These capillary jets are 
non-uniform in luminosity and variability within the half-opening angle 9j o f the outflows (distinct from 
the half-opening angle 9h on the black hole event horizon due to collimation, ILevinson & Eichlerl (| 1993b : 
van Putten & Levinsonl (120031 )). that reach their maximum at the interface with the environment or a sur- 
rounding bary on-rich disk or torus wind. Surfaces of maximal flux (9 ~ 9j) should hereby have maximal 
luminosity an d variability, the latter by the expected turbulent boundary layer between the jet and a colli- 
mating wind (Ivan Putten & Gup ta 2009), consistent with the observed positive correlation between the two 

dReichert et al.ll200lb . The duration of the capillary jets is determined by the lifetime of black hole spin. 

1/2 

Identified with the duration T90 of long GRBs, it gives rise to a spectral-energy correlation Ely oc E p T 9 q 
between EL and th e peak energy E p in gamma-rays in agreement with observations by HETE II and Swift 
(|van Puttenl l2008b) with a Pearson coefficient that is improved by including T90. 



The Swift event GRB060614 of duration T90 = 102 s introduces a new class of long GRBs with no 
supernova and with a late-time X-ray tail (XRT), which generally ex tends on a ti mescale of 1 - 1000 ks 
common to long (T90 > 2 s) and short GRBs (T90 < 2 s) events alike (Zhang 2007b. This suggests there is 
no apparent memory in the GRB remnant to the initial state giving rise to the prompt GRB emissions. The 
common endpoint of long and short GRB from rotati ng black holes a re black holes with slow spin with no 
memory to initial spin. Subject to late time accretion dRosswogll2007r) . it would produce very similar XRTs 
to long and short events. 
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3. Spindown of black holes and proto-neutron stars 

For black hole inner engines, we note that the spin energy of a Kerr black hole behaves remarkably 
similar to that of a spinning top in view of a ratio j^j = \ cos~ 2 (A/4) = 0.5 — 0.5858, where sin A = 
a/M = J/M 2 (0 < |A| < vr/2), that remains remarkably close to the Newtonian limit of 1/2. Surrounding 
by high density matter, most of the black hole luminosity Lh provided by E sp i n is expected to be incident 
onto the inner face of a torus by equivalence in poloidal topology to the magnetos phere of neutron stars, 
whereby its mass M and spin evolves in response to a torque T evolve according to (Ivan Putter! 1999) 



-M, T 



-J 



(1) 



as a function of the black hole angular momentum J and mass M. This interaction can also be seen by 
inspection of Faraday's eq uation for the elect romagnetic field (E, B) as a consequence of frame dragging in 
Faraday's equations (e.g. iThorne et al.1 (|1986)) 



V X «E = -8 t B + 4tT J m , 

where a denotes the gravitational redshift and Vj denotes the 3-covariant derivative, and 



(2) 



(3) 



in terms of a Lie-derivative of the magnetic field with respect to uj. The role of © in an inner torus 
magnetosphere around a black hole is schematically illustrated in Fig. Q] 

A derivation of J m based on the two Killing vectors of the Kerr metric is as follows. Let F = u A e + 
*u A h be the four-vector representation of the electromagnetic fie ld F ab associ ated with u b , u c u r . = —1, of 
zero-angular mome ntum observers (ZAMOs) ( Lichnerowitz 1967 ). Following ( Shapiro & Teukolsky 19831 : 
Thorne et al.ll 19861) . u = — adt and u = a _1 (k + cjm) in the Killing vectors k b = (dt) b and m b = {d^) b 
of the Kerr metric. Consequently, V c u c = 0. ZAMOs measure e b = u c F ac and h b = u c * F cb , where 

2 - -1). To evaluate V a * F ab 



0, we note 



e = (0, E 1 ) and h = (0, B l ), where * is the Hodge dual (V 

V a (u a h b - u b h a ) = C u h b + (V c u c )h b - (V c h c )u b , (4) 
(u c V c )h b — (h c V c )u b . Projected onto surfaces of constant coordinate time t (orthogonal to 



where C u h b 



(C u h) ± = a- 1 (d t B + C u ,B) 



(5) 



in the frame of ZAMOs, where L w i is with respect to uj % = urn 1 (where m l is not a unit three-vector). 
Next, we decompose V a = D a — u a (u c V c ), and note *(u A \i) a bcd = £abcdU c e d , and the acceleration 
(u c V c )ub = a _1 V(,a. Consider 



V"(e abcd n c e d ) = e abcd (D b u c )e d - e abcd u b a c e d + e abcd u c V b e d 



(6) 



The projection of the right hand side onto surfaces (r, 9, 4>) normal to u b satisfies en,cd{E> b u c )e d + lij^a?e k + 



e ibcd (D b u c )e d + a l e ijk V : > (ae k ), where e aijk u a = e ijk = VhA ijk with 



ayh over 
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the 3-volume \fh of spacelike coordinates using \jk, A123 = 1. Here, the first term on the right hand 
side vanishes, since D^u c is spacelike: u (Di,u c ) = by construction and u c DhU r . = by normalization 



u = — 1. The result is an extra current C u iB in Faraday's equation. See lThorne et all (119861) and references 



therein for an alternative derivation. 





Fig. 1 . — The interaction of frame dragging with magnetic fields is described by an additional Lie derivative 
Jm © in Faraday's equations (|2]). It induces poloidal current loops between rotating black holes and sur- 
rounding matter via an inner torus magnetosphere, here shown in poloidal cross-section with no-slip bound- 
ary conditions on the inner face of a torus and radiative (slip) boundary conditions on the event horizon of the 
black hole. The result is a potentially powerful spin-connection between the black hole (non-rotating, left; 
rotating, right) and the torus, wherein the black hole event horizon serves as a passive load. The induced 
poloidal currents mediate energy and angular mome ntum transfer by Maxwell stresse s on the event ho rizon 
(|Ruffini & Wilsonl 1 19751 : blandford & Znajeklll977l ) and the inner face of the torus (|van Puttenll 1999m . If 
the black hole spins rapidly, Lorentz forces due to the induced poloidal currents cause the inner face of the 
torus to spin up, equivalent to those in the magnetosphere of a spinning neutron star by which it spins down 
( van Putten & Levinsonl l2003 ) . 



Applied to an inner torus magnetosphere (Ivan Puttenll 19990 as in Fig. [T] we have 



07T 



z 2 sin 2 A 



(7) 



(z 2 + sin 2 A) 3 

where the inequality refers to an ingoing poloidal magnetic field with the orientation and sign of J m as indi- 
cated. In this process, the event horizon serves as a passive load while torus, subject to p owerful competin g 
torques acting on the inner and the outer face of the torus, develops forced turbulence (Ivan Puttenl l 1999). 
The complete evolution is now described by a system of two ordinary differential equations for the black 
hole lu minosity Lu and the torque Th by the ind uced Lorentz forces exerted (mostly) onto the surrounding 
matter (Ivan Puttenll2008al : lvan Putten et allboi lb 



l h = —M = k (n H - n T ) n T , t h = - j = « (n H - ^t) 



(8) 



parameterized by a spin down coefficient k given by the variance in poloidal magnetic field energy, the 
angular velocity of the torus and Qh of the black hole in their lowest energy state. Here, the angular 
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velocity of the inner torus magnetosphere is equal to Ut by no-slip boundary conditions on the torus and 
slip boundary conditions on the event horizon of the black hole. This output (Lh,Th) will be re-radiated 
in various emission channels by the inner disk or torus in a catalytic process, mostly so in magnetic winds, 
MeV-neutrinos and gravitational waves. 

In modeling long GRBs from rapidly rotating black holes , the open magnet ic flux tube is supported by 
an equilibrium magnetic moment of the black hole ( Waldfl974 : van PuttenlboOlh whenever the black hole is 
exposed to an external poloidal magnetic field. We propose a positive correlation between 6h and the radius 
of the torus described by a jet luminosity (Ivan Putten & Guptall2009h 



Lj oc QjjZ £j 



(9) 



= zM for a black hole of mass M, oc {Q,tRt) 2 ^{z) is the kinetic energy in the torus, 



where Rt 

and e(z) = -y/l — 2/3/z denotes the specific energy of matter in an orbit with angular velocity Qt at 
the ISCO (|Bardeenlll970h . With n = 2, Lj scales with the surface area enclosed by the torus, while the 



scaling to £k follows from a s tability limit on the poloidal magnetic flux energy that the torus can support 
(Ivan Putten & Levinsonl 120031) . The limited half-opening angle of 6h — 0.15 rad from the event horizon 
accounts f or the fraction of about 0.1% in true energies in gamma-rays relative to the spin energy of the 
black hole (Ivan Putten & Levinsonl l2003). 



For PNS, we recall that the e volution of their angu l ar vel ocity with spin aligned magnetic field derives 
from a spindown luminosity (e.g. Kalapotharakos et all (|201 lb ) 



L 



PNS 



oc n 4 



(10) 



as a function of its angular velocity Q. Here, we assume that gravitational-wave emissions are negligi 
ble, which othe r wise might derive, with considerable uncertainty, f rom a variety of dynamical processes 
(e.g. iRees et all Jl974h:bwen et all Jl998h : lcutler & Thornel (boO^ i lHowell et all (boO^ i Undersson et al. 



(|201 lh : iRegimbaul (|201 lj) ). The rotational energy, luminosity and spindown time for an initially maximally 
rotating PNS are 



E rot ~E Cl Le-lO 51 ^ 



M 



lAhM P 



ergs" 1 , r c ~^ = 30Bi- 6 2 s 



(11) 



for a magnetic field strength B = i? 16 10 16 erg, where 



E c = 3 x 10 52 erg 



(12) 



where i? c = 3 x 10 52 erg denotes the maximal rotational energy of a PNS. The luminosity Lpj^s is radiated 
off in a magnetic wind that may serve to drive an explosion and especially so following core-co llapse of 
relatively heavy progenitor stars, while the possibility of it driving a GRB (IMetzger et al.ll201ll) depends 



crucially on baryon loading that may be prohibitive during the initial hot phase of the PNS (IDessart et al. 
200il2008h . 
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4. Template light curves for spindown of black holes and proto-neutron stars 

The light curves of long GRBs show considerable diversity in durations, count rates, and variability, 
wherein uncertainties in photon count rates are relatively small given the large count rates that are often 
observed. Our focus is on extracting information on the secular (slow time) evolution of the inner en- 
gine as described in the previous section in the sense of an ensemble average. We do not intent to model 
detailed behavior in any individual burst as it may arise from, e.g., intermittencies in any of the (magneto- 
)hydrodynamic and radiation processes. Correspondingly, we set out to normalize the light curves of long 
GRBs to unit duration and integ rated photon counts, and average these to extract a normalized light curve 



(nLC) (Ivan Putten & Gu pta 2009). 



In the present study, we consider light curve model templates for GRBs from spindown of black holes 
by © upon solving numerically for ([8]) with the closure relations for given by 

Model A: Q T = Q ISCO , (13) 
Model B : n T = -Sl H , (14) 

Model A considers the torus to be at the ISCO, where Qrsno = ± 3/2 1 . for co-rotating (+) and counter- 



rotating (-) orbits, where z = r/M (Sh apiro & Teukolskylll983h . For an initially maximally rotating black 



hole, the ISCO (+) expands in size from M to 6M in the Schwarzschild limit of zero rotation for co-rotating 
orbits. Model B considers the torus to be further out. The templates are calculated by numerical integration 
of (H]) with closures (fl4l) with an initially maximal spin. The minor energy output in gamma-rays ([9]) are 
modeled subsequently as a function of the calculated evolution of the black hole. 

On the basis of ©,([81) and (fT4l . we calculate a model light curve Lj(t) of the minor output in hig h 



energy emissions, to represent the intrinsic light curve in gamma-rays detailed in lvan Putten & Gup ta (2009). 
Fig. |2] shows the result for closure A, where Lj(t) starts at a finite value Lj(to) > at the time of onset 
to for an initially maximally spinning black hole, and gradually increases to a maximum when a/M = 
0.8388 before decaying to a finite value as Qjj approaches Qisco- For an initially extremal black hole, the 
maximum is attained with a delay 

^16% (15) 

-^90/ A 

relative to Tgo of the model burst with Lj (t) — Lj (to) > 0, where max Lj (t) / Lj (to) = 2.27. With closure B, 
Lj(t) starts promptly at near-maximum, and rapidly decays with black hole spin. For an initially maximally 
spinning black hole, numerical integration of ([8]) shows an overall efficiency of close to 60% (equal to when 
a/M = 0.8 initially) for A and an overall efficiency of 35% for B, as follows by direct integration of © 
subject to the two alternatives (fT4l . These results point to dissipation of a major f raction of black hole spin 



energy "unseen" in the event horizon, creating astronomical amounts of entropy (|B ekensteinl 1 1 9731) . Black 
hole spin down by ((8]) is largely viscous. 

As shown in Fig. |2j the luminosity in template A is initially anti-correlated to black hole spin while 
a/M ~ 0.8388, associated with the increase of Qn due to the expansion of the ISCO. This anti-correlation 
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is consistent with a rec ent observation of no or a weakly negative correlation between jet luminosity and 
spin rates that are high ( Fender et al.ll201Ch . 




20 40 

time [s] 
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(e H /max(e H , 



Fig. 2. — Shown is the template A for the light curve in gamma-rays for an initially extremal black hole and 
the associated horizon half-opening angle Oh following ©, representing spindown of a black hole against 
matter at the ISCO. The luminosity initially rises to maximal luminosity as Oh increases with the expansion 
of the ISCO. At high spin rates 0.8388 < a/M < 1, therefore, there exists an anti-correlation of luminosity 
with black hole spin. Template A closely matches relatively smooth light curves such as those of GRB 
970508 and GRB 980425 (Fig. 1 in Ivan Putten & Guptal (12009k 



A light curve template for GRBs from PNS may be derived from their spindown luminosity (TTOb . 
Characteristic values for their rotational energy, luminosity and spindown time are those of (fTTTt . We thus 
have 



Model C: LAt) 



(l+t/T c 



(16) 



upon neglecting any gravitational wave emissions as mentioned above. 



In comparing model light curves A-C above with the light curves in the BATSE 4B catalogue in the 
time domain, we shall assume that photon count rate integrated across the four BATSE energy bands is 
proportional to the luminosity Lj(t). This assumption refers to the temporal behavior in any given burst, not 
across different bursts that further involves a diversity in parameters, e.g., baryon loading, that may affect 
the efficiency in prompt GRB emission processes. The template for Model A is shown in Fig. [2 



For the purpose of searching for evidence of spindown of a black hole or neutron star in long GRBs, 
we set out to consider long GRBs on average. To do so, we shall exploit BATSE catalogue of long GRBs. 



5. Matched filtering applied to the BATSE catalogue of 1491 long duration bursts 

BATSE produced by far the largest catalogue of GRB light curves. Since there is no report of any 
degradation of the BATSE instrument with time, this method ensures a completely unbiased and blind se- 
lection of events. It reveals a bi-modal distribution of short and long bursts, represented by Tqq < 2 s and 
Tqq > 20 s, respectively, as summarized in Fig. [4] 
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Fig. 3.— (Top.) Overview of the durations T 90 and T 50 in the 1491 long GRBs in the BATSE Catalogue, 
here sorted by Tqq. (Bottom.) The distribution is bi-modal in both durations, while the log 10 T50/T90 ratio 
(with a mean of 0.407) is subject to considerable scatter due to variability in the GRB light curves. 



The 1491 light curves of long GRBs comprise 531 and 960 events with durations 2 s < T90 < 20 s and 
T90 > 20 s, respectively. For the purpose of our study, we first consider the normalized light curves (nkj) 
of each burst (i = 1, 2, • • • N) individually as a function of normalized time r = t/Tgo in the interval 

/=[-l,3], (17) 

subsequent to filtering out high frequency fluctuations by applying a moving averaging over 40 samples over 
64 ms (corresponding to 2.56 seconds) to the BATSE light curve data. For the BATSE catalogue of 1491 
light curves with Tqq > 20 s, the results are summarized in Figs. H] 
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Fig. 4. — Compilation of the complete BATSE catalogue of 1491 light curves of long bursts sorted by 2 s 
< Tgo <1307 s. Each light curve shown represents the sum of the photon count rate in all 4 BATSE energy 
channels, is smoothed with a time scale of 2.56 seconds and is plotted as a function of time normalized to 
T90 in the interval of [—1, 3]. We note an apparent trend from relatively similar to increasingly diverse light 
curves with Tgo- 



We next apply matched filtering to each nlcj (i = 1, 2, • • • N) of the BATSE light curves in Fig. [4] 
It consists of a fit seeking maximal correlation between template and the nlcj by optimizing Pearson co- 
efficients, in which mean values of photon count rates are ignored. It comprises a three-parameter fit to a 
template light curve tLC following translation Ati and scaling at of time, as well as a scaling in photon 
count rate, for an optimal fit in the interval [—1,3]. This procedure may be seen to be equivalent to an 
optimal fit in the sense of least squares. Subsequently, the inverse of the transformation of the template thus 
obtained is applied to the GRB light curve, i.e., a shift —Ati to the origin and normalization a^ 1 in duration 
and /3- in count rate. A sample of the resulting nlcj is shown in Fig. [5j Each match to the template fLC(ifc) 
may be quantified by the standard deviation Uj of the residuals 



Pi (t k ) = nlCi(t fc ) - tLC(tjfe) (i = 1, 2, • • • ,N, k = 1, 2, • • • , v, t k e[-l, 3]) (18) 



where v = 1600 in our numerical evaluation. For Template A, the results are shown in Fig. [5] 
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Fig. 5. — (Top). Shown is a sample of the nlcj (blue) of the first 32 light curves from Fig. |4]by matching to 
a template (red) using the linear transformations of shift and scaling of the time axis and photon count rates. 
(Bottom.) Shown is the standard deviation (STD) of the deviation between the nlcj (i = 1, 2, • • • , 1491) 
and Template A and its probability density function. The STD increases with duration, resulting from an 
increasing variability with duration T90 with a shallow minimum around T90 = 10 s. 

An a priori estimate of the resolving power of our matched filtering procedure may be derived from the 
Standard Error in the Mean (SEM). In considering sub-samples of size N, such as in moving averages to be 
considered below, the tandard deviations of the residuals (1211 shown in Fig. [51 are 




We next set out to calculate the nLC obtained from various model templates. 



6. Extracting a normalized light curve (nLC) of long GRBs 

Since the nlcj are normalized, they are bounded on [-1,3]. The nlc, can now be averaged to obtain a 
normalized light curve (nLC) representative of all light curves of long GRBs in the BATSE catalogue: 

nLC(i fe ) = Isj^nlcifo) (20) 
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for the ensemble average of the nlcj(^) with standard deviation a(tk) at each discrete time t^ = — 1 + 
Ak/ue[—l, 3] (k = 1, 2, • • • u). For a final comparison between the nLC and the template, the nLC is given 
the same mean value in normalized count rate as the template. 



2s<T„„<20s 






-10 12 3 



— deviation 

- - ±3cj 



normalized time 



normalized time 





normalized time 



normalized time 



Fig. 6. — Shown are the nLC (circles) generated by model templates A-C (lines) for the ensemble of 531 
long duration bursts with 2 s < Tgo < 20 s (left) and the ensemble of 960 long bursts with Tgo >20 s 
(right) and the associated deviations for templates A-C. Here, the standard deviation a is calculated from 
the square root of the variance of the photon count rates in the ensemble of individually normalized light 
curves as a function of normalized time. 



In this process, baseline levels in background gamma-ray count rates are suppressed by matching tem- 
plates to data up to an arbitrary offset in count rate. By design, our procedure filters out short and intermedi- 
ate timescale fluctuations or modulations, that may result from intermittencies in accretion, the torus, shocks 
and turbulence in the outflow, leaving only the persistent evolution on a secular time scale of evolution of the 
black hole in each event. By the central limit theorem, (1201) hereby has a well-defined limit as N approaches 
infinity with asymptotically Gaussian behavior in fluctuations in the process of convergence. 

The nLC generated by the templates tLC of models A-C are shown in the left columns of Figs. [6l To 
quantify the uncertainty in the nLC resulting from the large but finite sample average, Table 1 lists the mean 
standard deviation a(tk) over tf. = [—1,3] in the subsamples of 531 individually normalized light curves 
nlcj(tfc) of the two groups 2 s < Tgo < 20 s and Tgo > 20 s. It expresses the degree of convergence in (1201) 
point wise, for each t^. The results in Fig. |7]show strong convergence in the nLC(tfc) at each t^ in the sense 
of a small standard error in the mean SE^ ~ 0.5%, which is less than the residual deviation S between the 
nLC and the template at hand as listed in Table I. The resulting Xred in (|22l) hereby effectively reduces to 
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Fig. 7. — Shown is the degree of convergence in the nLC expressed by the standard error in the mean SEj 
in (l20l with i =A-C. The three numbers refer to the interval / = [—1, 3], the interval / = [—1, 0] capturing 
the rise and the interval I = [0, 3] capturing the slow decay post-maximum. Because of the large size of the 
BATSE catalogue, the results are on average well below 1%. 



the ratio ~ <5/SE,;, measuring a residual error in the templates relative to the uncertainty in the observation. 
According to Table I, the residual errors are a few percent for Template A, which is more than adequate to 
distinguish it from various model alternatives. 



Table 1 lists the deviation between the nLC and its generating template, 

5(tfc)=nLC(t fc )-tLC(t fc ). (21) 

expressed in S[J] = STD ({S(tk)} tke j S j over the intervals / = [—1,3], / = [—1,0] during the rise and 
/ = [0,3] during the fall post-maximum in the nLC, and the goodness-of-fit in the match between the nLC 
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TABLE I. Summary on the nLC over tf.e[— 1, 3] (k = 1, 2, • • • , v = 1600) of two subsamples of lonj 
GRBs to templates derived from spindown of black holes (fT4l and PNS ( fT6l ) shown in Fig. [6] 



Quantity 



mean duration 
ensemble size 



Symbol 2 s < T 90 < 20 s T 90 > 20 s 



T 90 [s] 

N 



9.7 
531 



74.2 
960 



(sub-)interval 



[-1,3], [-1,0], [0,3] [-1,3], [-1,0], [0,3] 



convergences SE A [%] 0.455,0.461,454 0.619,0.731,0.587 

SE B [%] 0.474,0.507,0.463 0.0.581,0.654,0.561 
SE C [%] 0.405,0.393,0.407 0.469,0.466,0.470 



deviation 5 



Sa[%] 
Sb[%] 
Sc[%] 



goodness-of-fit x 2 red x\[%] 

xl[%] 
xl[%] 



4.0, 5.6, 2.9 
7.3, 12,4.2 
10.4, 10.0, 9.3 

69, 13, 35 
252, 860, 48 
1278, 1222, 1050 



3.2, 5.0, 1.8 
5.1, 8.8, 1.7 
6.6, 7.0, 4.7 

27, 65, 8.4 
93, 295, 10.7 
183,278,54 



to the template tLC by 

"nLC(t fc ) - tLC(ifc) 



N 

Xred[I\ = \ jJ| S <fc^ 



cr(t k ) 

where a(tk) is defined by the SEj(^) shown in Fig. f7]). 



(22) 



The results unambiguously show three trends in x 2 re fr pointing towards Template A, a prolonged decay 
in the nLC post-maximum on the interval I = [0, 3] and the sub-group of events with durations Tg >20 
s. Here, Template A is favored with a residual deviation 5a = 1.8%, whose x^[0, 3] = 8.4 is significantly 
smaller than the x^[0, 3] of Templates B and C. The residual deviation in Template A of a few percent is 
well beyond the intended accuracy of our model light curves. (More on this in the subsequent section.) 

The exceptionally large size of the BATSE catalogue allows us to elaborate in some detail on x 2 re d as a 
function of duration Tqq. Fig. [8]shows \ 2 f° r a moving average of ensembles of 60 light curves as a function 
of their mean T90. The results for template A show a clear break in convergence and in Xred at ^90 — 20 
s, or about 10 s in the comoving frame of the events based on their general correlation to the cosmic star 
formation rate. Quite generally, a number of processes may be present with time constants up to seconds 
in all long bursts, that are independent of the durations. Imprints of these processes in the nLC diminishes 
with increasing mean T90 of the ensemble at hand. The na confidence interval for x 2 re d associated with our 
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Fig. 8. — (Left.) Shown are the X 2 Te d\~ 1> ^] f° r tne °f Templates A-C as a function of mean durations 
T 90 of ensembles of 60 light curves with mean durations in the range Tg > 2 s. There is a noticeable 
transition in the match for Template A around Tgo ~ 20 s, which corresponds to a de-redshifted duration 
of about 10 s. The mean value of the Xred over 20 s < Tgo < 200 s are given in the legend. Overall, the 
discrepancy between the nLC and its generating template is smallest for model A. The largest contribution 
to Xredl~ 1> 3] stems from the onset in the interval [—1, 0]. (The dashed lines refer to 2a, 3a and 4a intervals 
of confidence.) (Right.) Shown are the x 2 [0, 3] for the prolonged phase of spindown post-maximum in the 
nLC in the interval [0, 3] of normalized time, which is parameter free for models A and C. The mean values 
1.43 and 1.74 of the xj! e£ j0, 3] for Models A and B, respectively, represent a 2.35<r and 4.05<r deviation from 
1. Spindown of a PNS (Model C) is essentially ruled out by Xred — ^-32 for all Tgo (a 12.7<t departure from 
1). (The dashed lines refer to 2a, 3a and 4a intervals of confidence.) 



moving average in time Tgo over subsamples of nlcj fe (ik = 1, 2, • • • 60) are (e.g. lAndrae et al.1 (120101) ) 



0.6349 < Xred < 1-3651 (n = 2), ( 23 ) 
0.4487 < Xred < 1-5477 (n = 3), (24) 
0.2693 < Xred < 1-7303 (n = 4), (25) 

represented by the dashed lines in Fig. [8] Fig. [8] further zooms in on the prolonged spindown phase in the 
nLC, post-maximum in the normalized time interval I = [0, 3] as a function of Tgo- Here, x 2 re d °1 Model A 
is within the 3<r confidence interval, while Model B is outside 4a and Model C is outside the 12cr confidence 
level. 
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7. Excess emission in gravitational-waves following a Hopf bifurcation 



The relatively small deviation between Template A and nLC^ of a few percent points to a sustained 
presence of an inner disk or torus that extends down to the ISCO for the duration of rapid spin of the 
black hole. In the presence of prodigious input from the black hole, this state requires cooling beyond that 
facilitated by magnetic winds and MeV neutrinos alone, as the following shows. 

An upper bound L* w on the luminosity in disk winds can be derived in the idealized limit in which 
all angular momentum transport Tjj in (O is radiated away to infinity by magnetic winds with vanishing 
emissions in MeV neutrinos and gravitational waves. We can estimate L* w by taking a weighted average of 
the local luminosity dL w = Qj,dA? associated with an electromagnetic vecto r potential A a supporting a 
net poloidal flux <I> = 2irdA ( k in a strip 2irrdr (cf. Ivan Putten & Ostriketi (120011) ). 



T* 



'' J?dE k ■ 



(26) 



assuming a constant ratio dA^/dE^ of poloidal magnetic field energy-to-kinetic energy in the disk. Here, 
the disk extends from z = rjsco/M to infinity. The associated energies follow by integration, Eh = 
Jq 00 Lffdt, E^ = f °° L* w dt following a prescription for the distribu tion dE k . A leading order approximation 
to dE k = 2irpHrdr is defined by the Shakura-Sunyaev solution (IShakura & Sunyaevlll973l) for a disk with 
radial density distribution, p, scale height, H, and perturbed angular velocity, SI/), satisfying 



r 

Po [ — 

Jo 



H = H n [ — 



n 



D 



(27) 



where the po, Hq and Qt denote respective scale factors. Resulting from a positive torque acting on the 
inner face provided by the central rapidly rotating black hole, we here have q > |. In this event, we have 

■kVL^pqHqtI 



E k = f 

J rn 



dE k = 7rfl T p H r 



s: 



- 2 *dl 



2q 



13 
4 



(28) 



provided that q > ^ (super-Keplerian motion). In this parameter range, dE k is a normalizable weight 
for calculating d26l ) that includes the condition q > \/3 (> ^) for a suffici ently slender inner disk or 
torus to be unstable to the n on-axisymmetric Papaloizou-Pringle instabilities (IPapaloizou & Pringldll984l : 
van Putten & Levinsonl l2003). 



In a core-collapse event, the explosion energy £$n produced by the momentum of relativistic magnetic 
winds and any conversion thereof into electromagnetic radiation onto the remnant stellar envelope from 
within satisfies £$n — ^PE W by momentum cons ervation, where den otes the velocity of the ejecta 
relative to the velocity of light. As a result, we have (|van Putten et al.ll2011bh 



E gw + E v = Eh — E w > Eh — E^ = 0' 



£sn, 



(29) 



where Eh / E^ is obtained by numerical integration shown in Fig. |9j relevant to the case of initial dimen- 
sionless spin values a/M > 0.4433 for which Eh — E^ is positive. 
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Fig. 9. — Shown is the ratio of the energy output of the black hole, Eh, to the maximum energy in disk 
winds, E* , in the process of spin down as a function of the initial dimensionless spin ratio a/M of the black 
hole. Beyond the critical value a/M ~ 0.4433, the excess energy Eh — E^ > is generally radiated 
in a combination of MeV neutrinos and gravitational waves. Observational constraints on the explosion 
energy in a relativistic core-collapse supernova attributed to disk winds hereby provide an upper bound on 
the gravitational-wave emissions. 



The excess energy Eh — E w in (|29l will be radiated away in an output E v in MeV neutrinos and 
and output E gw in gravitational waves. While the first is generic, the second may be seen to occur above a 
critical temperature. 



We recall that for infinitely slender tori, IPapaloizou & Pringlel (119841) derived an index of rotation q 
for the formation of non-axisymmetric instabilities in torus with radial angular velocity distribution Q(r) 
modeled by a power law, satisfying 



q > V3, fi(r) 



Ok I - 



(30) 



These instabilities open up angular momentum transport outwards by coupled wave-motion on the inner 
and outer face of the torus, and simultaneously s o for all non-axisymmetric modes to > 1. However, for 
astrophysical applications (e.g. IZurek & Go rski (1989)), we require an extension to wide tori beyond the 
singular limit of infinitely slenderness, described by finite ratios 0<5 = b/R<lof the minor radius b to 
the major radius R. Instabilities now occur for q > q c (6, to) > y3 for increasingly many m as the torus 
becomes more slender, recovering the bifurcation point q = v3 in the limit as 5 approaches zero. Fig. 2 
shows the neutral stability curves as a function of 5 for various to, a quadratic approximation to which for 
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\/3 < q c < 2 are Ivan Puttenl (120021) 
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showing that instabilities occur for q > \/3 for increasingly many m with decreasing slenderness 6. 

Heating produces thermal pressure (in addition to magnetic pressures that may be present), which in- 
creas es the rotation index, O-r(r) = {R/r) q ,£l 2 R = M/R 3 , about a point mass M, to (Ivan Putten & Levinson 



2003) 



1.5 + 0.15 



(— 

V 4M 



_6_ 

02 



T 



10, 



T10 =± 2L 



1/6 

v,52 



M 7 



0.1 Af, 







-1/6 



(32) 



associated with a ratio 5 = of the minor-to-major radius of the torus at temperatures T = TiolQ K. At 2 
MeV we note comparable Alfven and sound wave velocities in the presence of super strong magnetic fields 



(Ivan Putten & Gupta 



modes m = 1 , 2 are ( 



2009) . Critical values for the Papaloizou-Prin gle instabilities of the first two azimuthal 



van Putten 



lei 



120021 : Ivan Putten & Levinsonll2003h 
1 .73 + 0.27 j - , q c2 (5) = 1.73 + 0.27 



0.3260 



(33) 



where instability sets in for q > q c (6, m). The value of 8 in (l33l is assumed, and generally may depend on 
a variety of factors including the history of formation of the torus. As both thermal and magnetic pressures 
contribute to (l32l . it follows that a torus in suspended accretion is susceptible to pressure induced instabili - 
ties. Those that are non-axisymmetric strengthen by gravitational-radiation back reaction ( van Puttenll2002 ). 
whose luminosity is determined by the nonlinear saturation amplitude for unstable non-axisymmetric modes 
in the torus. An analytical estimate for the low-m multipole mass moments can be given for a flat infrared 
spectrum in MHD turbulence at the threshold of magnetic stability. The low-m mass inhomogen eities 
5m/ Mt will reach approximately the stability limit 6b /£& — 1/15 (Ivan Putten & Levinsonll2003r) . and 
produce broad line and quasi-periodic emissions in gravitational radiation. 

The onset of non-axisymmetric instabilities in th e torus for a critical energy input from the black hole 
(Fig. IT]) is reminiscent of a Hopf bifurcation (e.g. Keller ( 1986h ). here leading to a sustained emission in 
gravitational waves for the lifetime of black hole spin. We here identify it with the prolonged spindown 
identified in the nLC of the GRB light curves of §6. 



In ge neral, the observational consequences o f Models A and B follow from the equations of suspended 
accretion (van Putten & Ostriker 2001; van PuttenlboOlb . describing balance of energy and angular momen- 
tum flux from the black hole to the inner face of the torus (+) and that emitted to infinity (-) according 
to 



T+ = T_ + T GW + T„, ^+T + = Q_T_ + VlTTGW + Pu, 



(34) 
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□=0.07 Mpc 




Fig. 10. — (Left.) Diagram showing the neutral stability curves (solid lines) for the critical rotation index 
q c of buckling modes in a torus of incompressible fluid, which is an extension of the Papaloizou-Pringle 
instability to arbitrary ratios of minor-to-major radius b/R. Curves of q c are labeled with azimuthal quantum 
numbers m = 1, 2, • • •, where instability sets in above and stability sets in below. In q < 2, the m = mode 
is Rayleigh-stable. For q = 2, the torus is unstable for b/ R = 0.7385(m = l),0.3225(m = 2) and, 
asymptotically, for b/R = 0.56/m (m > 3). (Reprinted from Ivan Puttenl (120021) .) Here, the rotation index 
q as a function of slenderness b/R (dashed lines) is set by the temperature, here shown for T = TiolO 10 
K, Tio = 1,2,3 following (l32l . Heating is provided by the spin energy of the central black hole. The 



second window shows the critical MeV temperatures for which the index of rotation reaches the critical 
values q c (S, m) for the onset of non-axisymmetric instability. Generally, higher temperatures are needed 
for relatively wider tor i. This criterion may be relaxed by additional instabilities from magnetic pressures 
(|Bromberg et al.ll2008r) . The third window shows snapshots of the quasi-periodic gravitational wave signal 
during black hole spindown, scaled to a fiducial duration of about one m inute and plotted in frequency with 
a reduction by a factor of 1000 for illustrative purposes. (Reprinted from Ivan Putten ( 20083 ).) (Right.) The 
long duration negative chirp can be searched for using a dedicated time-sliced matched filtering procedure, 
here illustrated following signal injection in TAMA 300 noise dat a. It points to a s ensitiv ity distance of 
about 35 Mpc for the upcoming advanced detectors (Reprinted from Ivan Putten et all (1201 lh .^ 



where t + oc (Qh — ^+)> tgw an d r - denote the angular momentum fluxes in gravitational radiation and in 
magnetic winds, the latter with Lqw = Qtt~gw an d L w = £7_t_. Here, £l± denote the angular velocities 
of the inner and outer faces of the torus. In (l34l >. we have neglected energy and angular momentum loss 
matter outflow, that may be driven by the accompanying MeV neutrino flux. 

For model A, solutions to quasi-periodic emissions in gravitational-wave exist to d34l > whe n MHD 
turbulence has a flat infrared spectrum up to the first geometrical break m* ~ R/b (|van Puttenl (|200lh . 
neglecting r, y ), wherein luminosities scale w ith the energy in the magnetic field, and hence the kinetic energy 
in the torus (|van Putten & Levinsonll2003h . The gravitational wave luminosity is hereby determined self- 
consistently. At late times, as black hole spin down approaches the fixed point at wh ich its angular velocity 
approaches that of the ISCO in Model A, we have the gravitational wave frequency (Ivan Putten et al.ll201lb 



few = 595 - 704 Hz 



M 



(35) 



Vl0M Q 

where the frequency window of about 15% depends on the initial spin of the black hole. It predicts a fre- 
quency 1.5-2 k Hz at the end of binary coalescence of two neutron stars and considerably lower frequencies 
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in mergers of neutron stars with a black hole companion or core-collapse in massive stars ( van PuttenlboOSO . 



For model B, we consider (l34l > in the alternative limit with vanishing gravitational radiation. For slender tori 
described by small 5, [Q] = Q + — fl_ ~ qb/R = 2q5, the algebraic solution to (l34l ) with tqw = is 

ih 1 1 +0(8 2 ), pL = 26r 2 + 0(8 2 ), (36) 



n H 2i + ^(r 2 -i) v n l 



where T = 1/y 1 — v\j(?, assuming a symmetric partition of magnetic flux, in the inner and outer torus 
magnetosphere extending to the black hole and, respectively, to infinity and. We hereby arrive at Model B, 
representing a dominant emission in MeV neutrino emissions. 

As the preferred outcome of our matched filtering analysis, Template A points to a self-regulatory 
process of formation of multipole mass-moments: energy input from the black hole increases thermal and 
magnetic pressures which destabilize the inner disk or torus, while gravitational radiation provides c ooling 



which stabilizes by attenuation of the same. This mechanism is characteristic of a Hopf bifurcation (|Keller 



1986), describing the onset of stable oscillations at finite amplitude when a control parameter, here tem- 
perature, total pressure or rotation index, exceeds a critical value. The result is a stable state producing 
prolonged emission in gravitational radiation for the lifetime of rapid spin of the black hole. The luminosity 
in gravitational radiation scales with the energy in the poloidal magnetic field (the variance of the magnetic 
field strength), i.e., the m = component in the low-m spectrum of moments magnetic fields and mass 
of the inner disk or torus, that determines the spectrum in gravitational radiation. In the state of forced 
MHD turbulence in the inner disk or torus as provided by the input from the black hole, they are linked to 
dissipation at small scales that defines heating and radiation in MeV-neutrinos. In the state of forced MHD 
turbulence of the inner disk or torus, gravitational radiation and the MeV neutrino emission are inextricably 
linked representing, respectively, the infrared and UV-spectrum of its mass-motions. 

Very similar considerations may apply to spin down from rapidly rotating neutron stars except for the 
fact that its energy reservoir in angular momentum is about two orders of magnitude smaller than that of 
rotating stellar mass black holes expected to form in some of the most relativistic core-collapse supernovae. 
For this reason, these emissions are not considered here. 



8. Conclusions 

Core-collapse of massive stars produces black holes and PNS, where the apsherical relativistic events 
point to a rotationally powered inner engine. Some of these events are hyper-energetic, that can not all 
be attributed to spindown of a PNS. The association of some but not all long GRBs with core-collapse of 
massive stars points to inner engines with a hidden time scale of tens of seconds different from accretion. We 
therefore pursue a search for spindown of a black holes or PNS on the time scale of duration of long GRBs 
using three different model templates for the light curve of GRBs, defined by the spin of an initially rapidly 
spinning black hole by surrounding matter at the ISCO (A), further out (B) and spindown of a PNS (C). For 
each, we extracted a nLC from the BATSE catalogue of long GRBs by application of matched filtering. 
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By its large size, the BATSE catalogue allows the extraction of nLC's with an SE of less than one 
percent. This resolving power is adequate to distinguish various model alternatives considered in the present 
analysis. The match between uLCa extracted using Template A is within a few percent, whose prolonged 
decay post-maximum identifies the parameter free process of spindown in the Kerr metric. 

We find that nLC^shows a gradual rise in their light curve over about 16% of their duration according 
to (fT31) . that represents the initial phase of spindown of an extremal black hole (a = M) to a/M = 0.8388 
(Fig. A Xr e d analysis shows a best-fit for Template A relative the the alternatives B and C. As a function 
of duration, the x 2 re d g enera Uy improves with Tgo featuring a pronounced break at Tgo ~ 20 s in case of 
nLC^. Zooming in on the prolonged phase of spindown post-maximum in the nLC, the fit falls within the 
2. 35 c confidence interval, that should be contrasted with the x 2 Te d °f Templates B and C, which fall outside 
the 4a and, respectively, 12cr confidence intervals. Consequently, Template C, representing spindown of a 
PNS, is essentially ruled out. 

The break in the goodness-of-fit at durations Tgo ~ 20 s corresponding to a de-redshifted duration of 
about 10 s in view of the general correlation of long GRBs with the cosmic star formation rate. Beyond, 
it approximately levels off, providing evidence that their inner engines are normalizable. This is different 
from a standard energy reservoir on the b asis of a mean true energy output E~ ~ 10 51 in gamma-rays as 



originally envisioned by Frail et al.l (|200lh . in allowing for a diversity in the true energy in GRB-afterglow 
emissions. For instance, the brightest Swift events include GRB 050820A w ith a total energy output of 
Etot ~ 4.2 x 10 52 erg, i.e., £L and kinetic energies inferred from its afterglow (|Cenko et al.ll2010i) . 



The nLC refers to the shape in the normalized light curve of the inner engine in the limit of averaging 
a large number of individually normalized light curves. Inspection of in dividual BATSE light curv es reveals 
bursts with a smooth onset with remarkably good fits to template A (Ivan Putten & Gup ta 2009), as well 
as individual bursts with a prompt onset with remarkably good fits to either template B or C. Our nLC is 
a statistical result on the ensemble of BATSE bursts that is not intended to address features of specific, 
individual events or sub-classes of progenitors. Indeed, long GRBs may well derive from both core-collapse 
events and mergers and, if so, the onset of their respective light curves might be different, where the former 
is subject to the process of a break-out of jets t h rough the remnant st ellar envelope (associated with a time 
scale of about 10 s, e.g JBromberg et al.l (|201ir ): LNakar & S aril d20Tl| )), whereas the latter is not. 



A normalizable inner engine for long GRBs is consistent with having negligible memory on their 
formation history, especially when formed as near-extremal objects in core-collapse events. The distribution 
of their durations then depends largely on the poloidal magnetic field strength, that, for the first, further 
depends on the mass of the black hole and the surrounding disk. Thus, normalization of durations effectively 
collapses the intrinsic light curves to a unique light curve. Normalization will be effective when the durations 
are long relative to any of the small to intermediate time scales that are uncorrected to the duration of the 
bursts, such as time scales of accompanying (magneto-)hydrodynamic processes in the disk and the inner 
torus magnetosphere surrounding the black hole, in the evolution of the dynamo in a PNS, or time scales 
associated with the break-out of baryon-poor jets through a remnant stellar envelope. The former would 
include, for instance, dynamical time scales and time scales of instabilities, neither of which is expected 
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longer than a few seconds, while the latter is generally expected to be a few times the light crossing time of 
the progenitor star. 

The contemporaneous formation of ultra-relativistic capillary jets by frame dragging along the spin 
axis and spindown against surrounding matter are mediated by frame dragging. These processes combined 
are only beginning to be approached by numerical simulations. Earlier results show that ultra-relativistic 
j ets, sufficient to pro duce the GRBs that we see, can not be produced without general relativistic effects 



(Nagataki et al 



(INagataki et al 



2007). More recently, a positive correlation of jet formation with black hole spin is found 



201 lh . but a successful formation of ultra-relativistic jets awaits further developments. Our 



model A suggests that this may require pair-creation (which falls outside the realm of ideal MHD) along 
open magnetic flux tubes that extend from the black hole event horizon to infinity. A numerical simulation 
of a complete GRB light curve would further require taking into account black hole spindown against the 
surrounding matter, where the latter develops a state of forced turbulence, instabilities and multimessenger 
radiation processes. 

GRBs from rotating black holes give an attractive outlook for unification, as they appear naturally in 
various progenitor scenarios associated with core-collapse supernovae and mergers of neutron stars with 
another neutron star or companion black hole. Those involving rapidly rotating black holes can account 
for long GRBs with supernovae in CC events, and without supernov ae such as GRB 060614 in mergers 
events. Long GRBs are he rein spin powered, not accretion powered (|Woosleylll993l : IPopham et all 1 1 999: 



20081) 



Chen & Belaborodov 2007). The latter would force the black hole to continuously spin up (|Kumar et al 



, here in contradiction with the observed long duration decay in the nLC over about 84% of T90 post 
maximum. Mergers producing long GRBs with no supernovae would include bina ry coalescence of two 
neutron stars, as all of these produce rapidly rotating black holes ( Baiotti et al. 2008). Any merger scenario 
for long GRBs with no supernovae a priori rules out inner engines in the form of rapidly rotating neutron 
stars. 

Our identification of black hole spindown in the GRB light cu rves shows that m ost of the rotational 
energy is transferred to the surrounding matter, proposed earlier in Ivan Putted (|1999f ). which far exceeds 
that observed in the true energy output in the GRB -afterglow emissions. The result is a major output in 
MeV neutrinos and gravitational waves by catalytic conversion in the inner disk or torus, since this energy 
exceeds that which can be radiated off in magnetic winds. While the MeV neutrino bursts are difficult to 
detect even from sources in the local Universe, collectively they inevit ably contribute to the cosmological 
background in MeV neutrinos that may be observationally relevant (cf. INagataki et all (120031) ). 



In the process of black hole spindown, we identify a Hopf bifurcation giving rise to sustained non- 
axisymmetric instabilities in the surrounding matter for the li fetime of black hole s pin, in addition to any 
magnetic -pressure instabilities that may give rise to the same (|Bromberg et al. 2008). It lends credence to a 
quasi-periodic emission in gravitational waves for a duration of tens of seconds that may be detectable out to 
a distance of 35 Mpc b y the advanced ground based gravitational-wave detectors presently in development 
(Ivan Putten et al.ll201ll) . The anticipated output in gravitational radiation that may be dominant over MeV 
neutrino emissions, promising one step beyond SN1987A, whose > 10 MeV neutrino burst (compiled in 
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Burrows & Lattimen (119871 )) already exceeded the output in all electromagnetic radiation. It is therefore 



fortuitous that the goodness-of-fit of nLC^ exceeds that of nLCe- 
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